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Abstract 

The present article deals with Macroscopic Quantum Coherence resorting only to basic quantum 
mechanics. A square double weh is used to illustrate the Leggett-Caldeira oscillations. The effect 
of thermal-radiation on two- level systems is discussed to some extent. The concept of decoherence 
is introduced at an elementary level. Handles are deduced for the energy, temperature and time 
scales involved in Macroscopic Quantum Coherence. 
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I. INTRODUCTION 



Triggered by a seminal article^ written by A J Leggett in 1980, research into Macroscopic 
Quantum Coherence (MQC) has yielded impressive experimental,-"- theoretical^"- and even 
technological achievements^"—. The ideas developed in the last thirty so years by Leggett 
and his collaborators have not only changed the way we understand the relation between 
quantum and classical behaviours, but are also crucial in the future development of quan- 
tum computing. The present article aims at explaining the basic phenomenology of QMC 
resorting only to basic quantum mechanics. Thus, we believe this article can be of interest 
for any student who has attended at least a one-year course in quantum physics, and for 
faculty members committed to introducing students into contemporary research. 

In order to explain briefly what MQC is, let us consider a particle in a symmetric double 
well potential (SDWP). Figure 1 depicts an example of such a potential. In freshmen courses 
we have been told what to expect when the particle is in a high-lying energy level in a nice, 
analytical, potential such as this: for states for which the change in potential energy within a 
de Broglie wavelength is much smaller than the mean kinetic energy, the specifically quantum 
features of the behavior result negligible and the classical description becomes adequate.— In 
that sense, classical behaviour can be considered as a limiting case of quantum mechanics^. 
Suppose, nonetheless, the central barrier in the SDWP of Fig. 1 to be of macroscopic 
width. Then, the predictions of quantum mechanics and classical mechanics certainly clash 
for this system. A classical viewpoint would demand two distinct localized states of stable 
equilibrium, situated at —Xq and Xq, while quantum mechanics predicts an even probability 
distribution for the (non-degenerate^^) ground level state (which is, of course, the more 
stable stationary state.) In fact, the ground statefunction for such a potential would have 
to look something like Figure 2. 

Moreover, the (odd) eigenfunction of the first excited level (Figure 3), and indeed each 
one of the stationary solutions of a SDWP, necessarily has an even probability distribution. 

What Leggett predicted more than thirty wears ago, and what actually happens in exper- 
iments carried out in SDWPs of micrometric and nanometric typical lengths, is the appear- 
ance of a two-fold degenerate ground level E', with the system oscillating in an harmonic 
fashion between two eigenstates, \L) (Figure 4) and \R) (Figure 5) locahzed, respectively, 
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FIG. 1. An example of a SDWP potential, with characteristic double minima and central peak. 

at the left and right of the central barrier. At ground level, the position expectancy value 
oscillates in the accordance with: 

(x)o(t) = (x)o(0)cosa;i. (1) 

This phenomenon, the so called Leggett-Caldeira oscillations, is closely related with the Rabi 
oscillations of atomic physics. It is explained as the result of the purported ground level E' 
resolving into a true ground level 

E^^ E' -hw/2, (2) 

endowed with an even non- localized eigensolution |+), and a first excited level 

E_^E' + hjj/2, (3) 

endowed with an odd non- localized eigensolution |— ). When a quantum system tunnels 
periodically trough the barrier of a SDWP with a central barrier of macroscopic length, we 
have Macroscopic Quantum Coherence. 

The states \R) and \L) have, each one on its own, a definite value of a macroscopic 
property (namely, the property of being localized at the left or the right of the barrier). 
At the same time, \R) and \L) are linear combinations of the states 1+) and |— ), which 
cannot be said to be localized. In order to understand Leggett's original motivation, notice 
the analogy between macroscopic SDWPs and Schroedinger's cat: the celebrated pet can 
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FIG. 2. A rendering of what a ground-level eigenfunction (solid curve) should look like for a SDWP. 
The potential is shown as a dashed curve. 

be in any of two different macros copically distinguishable states (let us say, for a live cat 
and \l/o for a dead one) just as a particle in a SDWP. If any of these macroscopic systems 
obeys the laws of quantum mechanics, then it could be prepared in linear combinations 
that lack a sharp, well defined, value of the macroscopic property. Examples of these linear 
combinations are the |±) states of the SDWPs, and the "neither dead nor alive" states 



of the cat. Thus, a more general definition of MQC is simply: the quantum superposition 
of distinct macroscopic states. Long time before the year of 1980, macroscopic quantum 
phenomena had been discovered: superconductivity in 1911, and superfluidity in 1937. Yet 
it remained for Leggett to identify the conditions necessary for a quantum system to present 
macroscopically distinguishable states.- 

Some twenty years elapsed between Leggett 's proposal and a credible experimental 
confirmation^'^ of MQC. One of the main reasons for this delay lies in the fact that the 
phase coherence of the |±) states is rapidly lost due to the interaction of the system with 
its surroundings, so the system collapses into one of the localized states before one period 
of the Leggett-Caldeira oscillation is completed.-"^ 

MQC is relevant not only from the purely theoretical point of view. A physical qubit is 
a two-level system considered as a piece of hardware. And, as we shall see in the following 
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FIG. 3. Th first excited level eigenfunction (solid) of a SDWP (dashed). 



pages, at least some SDWPs can behave as effective two-level systems at sufficiently low tem- 
peratures. Quantum computing (an area with impressive software development, but little 
hardware to show) requires qubits to interact with one another without loss of coherence, for 
fairly long times, even at fairly high temperatures. Thus, the study of two-level dissipative 
systems, to which Leggett and collaborators made far reaching contributions when delving 
in the foundations of quantum physics, has revealed itself crucial for people in the vanguard 
of technological development .-i^ 

The rest of this article is structured as follows: in section HIl we discuss the spectra of a 
family of symmetric double square well potentials, and the conditions under which a member 
of this family can be considered as an effective two- level system. Next, the properties of two- 
states systems arising from SDWPs are discussed in section lllli We then go on to examine 
in IIVI how thermal radiation, by throwing the system into higher energy levels, renders 
the two-level model inapplicable. In section |V] decoherence is introduced in elementary 
terms, and its relation with dissipation is briefly discussed. Handles for the time, energy 
and temperature scales involved in MQC are derived from our toy model in section IVI[ 
Finally, conclusions are laid down in section I VI II 
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FIG. 4. The {x\L) state (shown sohd) locahzed at the left of the SDWP (shown dashed) barrier. 
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FIG. 5. The state (shown sohd) locahzed at the right of the SDWP (shown dashed) barrier. 

II. SYMMETRIC DOUBLE SQUARE WELLS 

Leggett resorted to quasi-classical considerations when stating his original proposal.- 
Also, the WKB approximation has been applied to double well potentials by Landau and 
Lifshitz,— and more recently, in this Journal, by others.— Here will take a different point of 
view, avoiding all together quasi-classical approximations, by considering a particular family 
of double infinite square well potentials as approximations to actual, analytic SDWPs. Our 
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procedure will later allow us to get some reference values on the energies, temperatures 
and times involved in MQC The following family of piece-wise-constant potentials will be 
considered: 

oo if X < —a — b, 
if — 6 > X > —a — b, 
Uh(x) = ( k if b>x>-b, (5) 
if6 + a>a;>6, 
oo if X > 6 + a. 

Potentials of this kind have previously been studied in a different context, and it has been 
showni^ that, if all other parameters held fixed, levels -E'2n+i and coalesce as /c — i- oo. 
Here we shall consider the barrier height A; > as a fixed number, although "big" in a 
sense that will be readily clarified. This, in order to keep the gap between the ground and 
first excited levels sufficiently small. We shall also take the width of each one of the lateral 
valleys, a > 0, as a fixed value unless otherwise stated, leaving free the only other parameter, 
that is the barrier half-width 6 > 0. 

One of the two main objectives of this Section is to obtain a global lower bound on the 
energy gap between the first and second excited levels in the Ui, potentials. Just as important 
to our ends, we will learn on this Section that there is a "running" upper bound (dependent 
on the value of b) on the gap between the ground and first excited levels. The consequences 
of this to facts, which are vital to the rest of the article, are explored in Sections III, IV and 
VI. 

To be sure, non of the Ub is continuous, yet they share the most prominent features of a 
SDWP, namely, they are even potentials with completely bounded, non-degenerate, spectra, 
as can be shown from boundary conditions. If instead of two minima, the Ub have two 
non-overlapping regions of minima {viz. {—a — b, —b) and (6, a + 6)), this distinction will 
prove to be quite unimportant. 

Also from boundary conditions (or from more abstract, symmetry considerations) it is 
readily seen that the levels in the spectrum of any of the Ub are classified according with 
parity, just like it happens for a continuous even potential: 

i'lnA-^) = ^-in.bix), 72 = 0,1,..., 6 G (0, CX)) (6) 

and 

^2n+i,b(-a;) = -V^2n+i,6(a;), n = 0, 1, . . . , 6e(0,oo). (7) 



—a — b 



-b b 



a + b 



FIG. 6. A typical member of the Ub{x) family of potentials 

Let us focus on the discretization conditions below the level of the central barrier {E < k). 
Prom the boundary conditions we get, for even states: 



r— V^mE^ r — J2m{k- E2n) 

- \/E2n cot a ^ yjh- E2n tanh b- , 

while odd levels below the barrier level have to comply with 

r- ^/2mE2„+l r ^1m{k-E2n+x) 

- V -C'2n+i cot a = \Jk- £/2n+i coth 6 



(8) 



(9) 



Notice how the first of these two conditions can be written in the form: 

g{E2n) = hb{E2n). (10) 

and the second can be rendered as: 

5'(-E2n+l) =i6(-S2n+l), (H) 

with the meaning of g,hh and jh being obvious from the context. 

Both of these two last equations are depicted in Figure 7, from which it can be seen that 
there exists an upper bound B, given by 

such that the ground and first excited states have to comply with 

J < Eo < El < B, (13) 



no matter the value of b. Obviously, there can be no levels below the barrier unless k > B/4. 
We shall only consider potentials for which the condition: 

k» B (14) 

is met, so that we will always have at least two levels below the barrier. Indeed, the number 
of levels below the barrier increases with increasing quotient k/B and, more importantly, 
as B is independent of k, condition ( I14p warrants that the gap between the first two level 
is always small. It is not difficult to generalize f lT^ starting from and and definition 
f lT2|) . The result is that: 

{n+^fB<E2n,k<E2n+i,k<{n + iyB, n = 0,1,2, ... N, (15) 

if the level 2N + 1 is still below the barrier. 
From inequality f|T5|) it follows that 

^2n+2 - ^2n+i > + 5/4)5, n = 0, 1, 2, . . . , (16) 

if level 2N + 1 is below the barrier. We then have that the gap between the ground and first 
excited levels will always be less than the gap between the first and second excited levels: 

E2- Ei>^B >^B > El- Eq. (17) 

But we can do much more better than that. Indeed, in Appendix A it is formally proven 
that for any given number 5 > there exist a value h{5) > such the gap between the 
ground and the first excited level of a Ub potential will be less than 6, that is 

Ei~Eo< 6, (18) 

if 6 > b{6). In other words, if we choose the barrier length big enough, then we can make Eq 
and El as proximate as we want, while there is a lower bound for the gap between E2 and 
El which is independent of the value of this length. This will allow us to find examples of Ub 
that will work as effective two-state systems for the lowest-lying energy levels, as illustrated 
in Figure 8. 

Finally, there is one more inequality that can be derived from ( ITSl) and that will prove 
useful in section IV. This inequality is: 

E2-Ei<^B . (19) 

Let us stress that relations (fT3|) . (fTTI) and (fT9l) are verified for each Ub regardless of the value 
of b. 
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FIG. 7. Graphical solutions of transcendental equations ([8]) and ([9]). Depicted, functions g{E) 
(solid), hb{E) (squares) and jb{E) (crosses). In all cases k = AOB and b = 0.8a. In this example 
only the ground level and the three first excited levels are below the barrier height k. 




FIG. 8. Doubling the barrier width produces a dramatic decrease in the gap between the ground 
and first excited levels. Shown: function for b = 0.8a (squares) and b = 0.4a(circles), and 
function for b = 0.8a (crosses) and b = 0.4a (asterisks). In all cases k = 40B. 



III. TWO-LEVEL SYSTEMS WITH REFLECTION SYMMETRY 



In the preceding section we have proven that there are Ui, potentials for which the gap 
between the ground and first excited energy levels is much more narrow than the one between 
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FIG. 9. Solid curve: the ip+{x) ground- level eigenfunction (even), obtained through computer 
assisted numerical analysis. Dashed: the corresponding Ub potential. In this example k = 145 
and 6 = 0. 

the first and second excited levels. Consequently, for low energy expectancy values, a particle 
in one of such potentials acts as an effective two-level systems.— 

In the rest of this section we shall consider a fixed Ub that behaves as a two-level system, 
and drop the b. 

Consider now the non- stationary solutions ipi and ipR that one obtains from the linear 
combinations 



and 



exp ( -^-^j i'oi^) + exp (-^-^ ] i'li^) 



exp ( -^-^ j i'oix) - exp (-^^ ] i^ii^) 



(20) 



(21) 



These states have no definite parity, but instead one is the specular image of the other: 



(22) 



as can be seen from equations (E]), ([7]), (120|) and (12T|) . 

The position expectancy value for this states is calculated from ([6]) in a straightforward 
manner: 

{x)L{t)= - (x)ij(t) = {i)o\x\ipi) COS 1, (23) 
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FIG. 10. Solid: the (odd) eigenfunction of the first excited level. Dashed: the Uh potential to 
which this solution corresponds. In this example k = lAB and b = 0.2a. 



as is the energy expectancy value: 



(H), = = :5l±£i. (24) 



Comparing ( l23l) with ([T]) and ( l24l) with ([2]) one may be tempted to make the identifications 

E' = {H)l and oo = El^El^ (25) 

from which (|3]) would follow, so that the states of f l20|) and fl2T]) could be interpreted as the 
localized states observed in the experiments, and V'o and ipi would correspond to the true 
ground level and the first excited state E^. That is, it would be cogent that 

(x|L) = ^i(x), {x\R) = tljR{x), = ^o(a;), = ^i(x). (26) 

In this interpretation, however, there is no room for transitions. Indeed, the complete 
Schroedinger equation for a U potential, which reads: 

~2^d^^'''^^^ t/(x)^(x, t) = th^ix, t), (27) 

predicts that if the system is initially prepared in the state ipiix^t = 0) at time t = 0, 
then it will remain in the ipLi^jt) state for t G [0, oo) (which is a sophisticated way to say: 
forever). This is just consequence of PDE's theory. 
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FIG. 11. The localized tpzixjt = 0) eigenfunction (solid) of the potential (dashed) for b = 0.4a 
and k = UB. 




FIG. 12. The localized il)R{x,t = 0) eigenfunction (solid) of the Ui, potential (dashed) for b = 0.4a 
and k = UB. 

Instead of periodical transitions between two different states, our equations predict the 
existence of a unique "oscillating" state, because iPr{x, t) is a time-displaced replica of 

i/jiix, t+-) ^ lexp ( - I—) iiR{x,t) , (28) 
13 



where fl stands for 

n = {Ei + E2)/2h (29) 
and u is as in (125!) . One arrives at this result directly from (!20|) and (12T]) after some algebra. 

A. Flip-flops and Leggett-Caldeira osciUations 

Let us start from what we know happens in actual experiments {i. e. the existence of 
an observable degenerate ground level) and proceed to deduce from there the perturbation 
needed to achieve such degeneracy. The SDWP Hamiltonian H is represented by the matrix 

H= " ] (30) 

in the symmetry-respecting basis formed by the eigenfunctions ipQ and ipi. Let us consider 
another Hamiltonian, if', represented by the matrix 

, ^ ( E' q\ 

e' = OH'O"^ = (31) 

in the symmetry-violating basis spanned by ipi and ipR. Here, O stands for the unitary 
operator which transforms ipo into ipi and ipi into ipR, thus: 

Notice that, as the rhs of (|3T1) is proportional to the identity matrix, then: 

, , 

H' = . (33) 



E' 



Thus, the perturbation is represented by: 



. , 

W = H' - H = . (34) 

-f)hj/2 

Now, in the basis spanned by ipi and ipR the things look quite different. Indeed, we have 
that: 

e = oeo"^ = ( ^ -?«^/2 

-nw/2 E' 
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and most importantly: 

I tvjjl2 , , ^ 

W = OWO"^ = . (36) 

\^ nuj/2 

So that the perturbation has no diagonal elements. This means that zeroth order corrections 
are strictly null for the perturbation, and, further more, that the off-diagonal elements are 
equal. 

To be very clear, let us write the eigen-equations for each one of this distinct systems. 
For H we have 

Hi;o{x,t) = ih^{x,t)=E^i,o{x,t), H^,{x,t) = th^{x,t) = EiM^^t), (37) 
while H' responds to 

H'ijL{x,t) = th^{x,t) = E'ijL{x,t), H'ijR{x,t) = zh^{x,t) = E'i:R{x,t). (38) 

Now, let us consider W as the initial, unperturbed, Hamiltonian matrix, and 

- W = -OWO^^ (39) 

as the perturbation, so that EI is the final, perturbed, Hamiltonian matrix. Then we can 
show that the ipL{x,t) and ■ipR{x,t) states transit from one another in Rabi style. Indeed, 
resorting to the time-dependent perturbation formalism,—!^ we write, for a general state 

^{x,t) of m, 

^{X, t) = CL{t)ijL{x, t) + CR{t)^R{x, t), (40) 

in order to obtain the equation 

which is equivalent to the 2x2 system of coupled linear equations: 

dcL hhj dcR hw . 

in—— = cr ,ih—— = cr. (42) 

dt 2 ^ ' dt 2 ^ ^ 

By uncoupling this system we get the harmonic oscillator equation 



df^ 4 
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CL (43) 



and a similar equation for cr, so that 

CLit) = sin {ut/2 + 0) CR{t) = cos {ut/2 + 0) , (44) 

where is a constant that can be elucidate from initial conditions. The probability of finding 
the particle in the state ipi is given, according to this last equations, by 

Pl (t) = sin' (a;t/2 + 0), (45) 

and the probability of finding the the particle in the ipR state is 

Pnit) = 1 - PL{t). (46) 

This is a particular instance of the Rabi oscillation, and this case is resonant due to the 
degeneracy of the "initial" Hamiltonian H'. But the "perturbed" Hamiltonian H is nothing 
else than the SDWP Hamiltonian of equation fl27|) . 

Now, equations and (jlSD predict the "flip-flop" between the stationary states iPr{x) 
and so that, if the system is initial prepared in the state 

ip{x,t = 0) = Mx), (47) 

then we will have a 100% certainty to find it in state iPr{x) at times t = 5? 77? tj • • ■ 
and a 100% certainty to find it in state '0'l(x) at times t = 775 775 77 • • •• And this last 
result is consistent with equation ( !28|) . Thus, we are in the presence of two different (yet 
not contradictory) descriptions of one and the same phenomenon: if H is considered an 
unperturbed Hamiltonian, with complete stationary solutions tpQ{x,t) and '0i(a;,t), then we 
have an "oscillating" non- stationary solution ipiix, t). If, on the other hand, H is considered 
to be the result of a perturbation acting on the degenerate Hamiltonian H\ we then get 
flip-fiops between the complete stationary solutions of if', that is: ipL{x,t) and ipR^x.t). 

In this manner, we obtain the periodic transitions (the zero point Leggett-Caldeira oscil- 
lations) observed in so many experiments. Notice that this transitions occur in the absence 
of external fields, thus without emission or absorption. 

IV. THERMAL RADIATION 

Due to the fact that no quantum system can be completely isolated from its environment, 
in any realistic description the Schroedinger equation must be supplemented with terms 
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that describe the interaction between the system and its surroundings. But there are very 
different ways to describe this interaction and its results, depending on the time and energy 
scales involved, and the complexity of the analysis. Here we shall discuss the absorption- 
induced transitions by which the system is thrown into high-lying energy levels, rendering 
the two-level model inapplicable. The main result from this discussion will be a limit on the 
temperature at which Caldeira-Leggett oscillations can be observed. 

A. Oscillations near resonance 

Now, oscillatory behavior is to be expected not only for the the resonant, exactly degen- 
erate, Hamiltonian matrix W . Indeed, it would not be realistic to expect Leggett-Caldeira 
oscillations only in perfectly isolated systems. Consider a harmonic perturbation of the 
SDWP matrix Hamiltonian EI of equation (130|) . that is, a perturbative term of the general 
form 




(48) 



and let us focus on the particularly simple case for which 




(49) 



so that the perturbative term can be written down as 




(50) 



We then again resort to the time- dependent perturbation formalism, and write 




(51) 



in order to obtain the equation 




(52) 



where 03, defined by 




(53) 
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represents the perturbation in the interaction picture, and in our particularly simple case 
reduces to 



/ 



QJ(t) = A 







exp {i{uj' — uj)t) 



\ 



v 



(54) 



exp {—i{uj' — uj)t) 
so that equation (!52|) is equivalent to the 2x2 system of coupled ODEs 



ih—— = A exp 
dt 



i(uj' — uj)t 
It can be checked by hand that 



Ci, in—— = A exp 
dt 



in' 



— — uj)t 



Co. 



co(t) = exp(i(]'t/2)<( cos{Rot) - ——sm{Rot) 

2Rq 



and 



Ci(t) 



-iRi 

Rn 



exp{-in't/2) sm{Rot) 



(55) 



(56) 



(57) 



provide a solution for the initial conditions co{t = 0) = 1, ci{t = 0) = 0. Here we have used 
the following shorthand 



Rn 



{A/hf 



Q! = oj' — oj and Ri = A/h, 



(58) 



which lead to what is known as Rabi's formula,— namely: 



Ri 



sm 



{Rot); 



Poit) = l-Pi(t). (60) 

It is not difficult to find the expressions for Piit) and Pnit) for this particular choice of 
A. We omit these, as they are not particularly illuminating. Let us just point out that in 
all instances Pr and Pl are oscillating functions of time, although they are generally not 
periodic. If one wishes to describe periodic Rabi oscillations in the R and L states, one 
should take, instead of (H9ll . 

A = - 

2 1^1 -1 

as the natural choice for A. By doing this one obtains expressions completely analogous to 
dSni) and dSO]) for Pl and Pr. 



(61) 



(62) 
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PR{t) = l~PL{t), (63) 
where the new frequency of the oscillation is now given by 

i?'o = ^J{A|hf+{uJ'|2f (64) 

The main conclusion of this subsection is thus, that the Leggett-Caldeira can survive the 
influence of an environment on the particle in a SDWP under certain circumstances. 



B. A limit on temperature 

An important result from perturbation theory tells us that for harmonic perturbations 
the time-dependent transition amplitude, c„_>m(i), between two given eigenstates of the 
complete Hamiltonian H is given by:— 

, , , 1 - expzf ^'""^" - uj')t , , , 1 - expzf ^'"-^" + u')t , ^ 

As a consequence we get that, if the system is to stay in the two lowest lying levels, then 
the perturbation must meet the condition: 

E2- El 

uj < . (66) 

h 

Otherwise, the perturbation would excite the system to higher levels with non-negligible 
probability. This gives a limit on the temperature at which the system behaves like a low- 
lying two-level system. Indeed, recalling Wien's law for blackbody radiation, we get that 
thermal radiation at a temperature T will have a maximal contribution of frequency uo' when 
condition 

= (67) 

is met. (Here, T stands for the temperature of the radiation, hw is Wien's constant, and c 
the velocity of light) Thus, if V(x, t) is somehow to represent thermal radiation, and if the 
perturbed Hamiltonian H' is to be described as a low-lying two-level system, then we must 
have: 

T < (68) 

In so many words: for each system there is a limit temperature above which the two- 
level system description is inapplicable, and zero-point Leggett-Caldeira oscillations become 
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overshadowed by other transitions. Moreover, from inequahties (fT7|) and (fT9|) we get : 

TB{a,m) < ^ "^^ . "^^ < 3TBia,m) (69) 
/vrc Ai 

with this global bound given by: 

The meaning of expressions (l69l) and (!70!) is the following: consider a family of double 
rectangular barriers, with a fixed m, a and A;, but free barrier width. When exposed to 
thermal radiation, there is a temperature Tg for the radiation above which the Leggett- 
Caldeira oscillations are overshadowed by other transitions in least some the systems, and 
at temperature 3Tb the Calderia-Legget oscillations are surpassed by other transitions in all 
of the systems. 



V. DECOHERENCE AND DISSIPATION 



We begin this section with a simplified exposition of mixed and pure states and the 
density matrix formalism as found in Landau and Lifshitz,— to move on next to an also 
simplified rendering of some of Leggett's original argumentation. After that, decoherence is 
defined, an its relation with dissipation is briefly discussed. 

The interaction of a system {&) with its surroundings (2) can be taken into account 
by considering a bigger isolated system (il) which encompasses both & and (E (that is: 
il = © [J (£). The state of this new, all including, system, il is described by a state function 
\E'(j, ^) that depends both on the coordinates of & (the j) and the the coordinates of its 
environment (the ^). The total Hamiltonian Ht acting on ii can always be written in the 
form: 

Ht = H + H^ + XHj (71) 

where H depends only on the j and their generalized momenta, depends only on the ^ and 
its momenta, and depends on both types of coordinates. We shall take the approximation, 
that H is the Hamiltonian of & when isolated, and that Hj alone models the interaction 
between & and 

In principle, there can happy instances in which "^{j,^) could be written as the product 
of two states functions: 

'^{j,o = mm (72) 
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but this does not need to be the case. States that can be written in the form (1721) are called 
pure states in the literature. States that are not pure are said to be mixed. 

In order to illustrate this let us consider the case in which both the original system and 
its surroundings can be represented as two-level systems. If the isolated Hamiltonian H has 
eigenfunctions and ?/'_: 

Hil^i = ^±^± (73) 
and if 0q, and 0/? are the eigenfunctions of Hg, i. e. 

He4>a = Ea(l)a , He(l)p = Ep(j)p, (74) 

then some examples of pure states are: 

^(^+0/3 + V^-0/3) = 75(^^+ + V^-)0/3, 1^-0/3 + ^^-0a = ^-{\(i>P + ^0a) 

and 

On the other hand, as instances of mixed states, we can provide the following: 

--^(^_0„ + , -^(0'_0^ + i)+(t)a) , and -^(0'+0„ + ^+0/3 + ^-^a)- (75) 

The density matrix formalism was developed to treat systems that (like ii) can present mixed 
states. The density matrix p allows us to calculate the expected value (/) of any observable 
f{x,Px) that depends only on the coordinates and momenta of &: 

if) = Tr(/p). (76) 
The elements of the density matrix p of a state \l/(j, ^) of are defined as: 

p,,, = S^^*ij,0^if,0, (77) 

where stands for the sum over the discrete ^ (if any) plus an integral over the continuous 
^ (if any). In the case of our 2 x 2- level system, expression (1781) reduces to 

PjJ' = + J> / = ± • (78) 

The diagonal elements of density matrix, of the form called populations, while the 

off-diagonal elements {i. e. the elements with j ^ j') are known as coherences. 
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Suppose now that the 50-50 hnear combinations 

i^L = ^ (^+ + , = (V'+ - V-) (79) 
are eigenf unctions of a macroscopic observable M, let us say: 

MijL,R = PL,Ri^L,R, (80) 

and take then the mixed state given by 

^ = CLi'L(pa + CRi)R(j)p. (81) 

The density matrix associated with (ISTj) is written as 

pH'^'V", I (82) 



\CR 



2 



in the {t/'l, V'i?} basis, as can be seen from (178|) . so that according to ( !76|) the expected value 
of any observable / pertaining to & yields the value 

{f) = \cL?fL+\cR?fR, (83) 

where fi and //? are the expected values of / in the pure states 

^L = ^L<^a and mji = iljR(t)p. (84) 



The point of this discussion is that the same result fl83|) is obtained if we make measurements 
on an ensemble of il systems all in state \E', or if the same measurements are made with an 
ensemble of it made up of a statistical mixture of the pure states ^ l and ^ r, in proportions 
|clP and \cr\'^. If it were to be held true that only ensembles of the type flHTj) could be 
prepared for il, then it could be argued that property M has a sharp value for each element 
of the ensemble, and that a measurement done on a particular element only removes our 
ignorance on its value for that particular system. Clearly, this opens the door for hidden 
variable theories. To put it succinctly: in this interpretation each one of the Schroedinger's 
cats in an ensemble of such felines would be either dead or alive, and never in superpositions 
composed of both dead and alive states. Only the behaviour of the ensemble would be 
quantal, its individual elements being essentially classical. 
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It is patent, on the other hand, that the pure state 

^+ = (85) 

cannot be written as a mixed state of the form (IHTj) and that its corresponding density matrix 
cannot be diagonal in the {L,R) basis, unhke ( l82l) . The impossibihty of the simultaneous 
diagonalization of the density matrices of all possible states of a system & is then a strong 
evidence of the true quantal behaviour of such system, as opposed to the behaviour required 
by hidden variable theories. Thus, for a system S to be classical in any sense of the word, 
the coherences, i. e. the off-diagonal elements, must be absent from the density matrix for 
each one of its possible states. This conclusion is generally valid, even if we resorted to the 
most trivial case in order to illustrate it.- 

Decoherence can be defined as the decay of the off-diagonal elements in the density matrix 
as a result of the interaction of the system with its environment. Therefore decoherence 
allows a system to behave as quantal when isolated and as classical when the coupling with 
its environment is "sufficiently effective." This is nowadays considered a plausible mechanism 
for the emergence of classical reality from a quantal substratum. 

In most practical applications, the environment (£ has a very large number of degrees of 
freedom (say of the order of the Avogadro number) and not just one, as in the example we 
have used. Thus 11 is usually a thermodynamic system, so that the full toolbox of quantum 
statistical mechanics needs to be marshalled in order to describe its behaviour. In this case 
the interaction between & and € (interaction known as quantum dissipation in this context) 
involves the relaxation of the thermodynamical variables of il towards thermal equilibrium, 
and not only decoherence. 

Various models have been proposed over the years for the environment (or bath) but 
one of first and most successful is the spin-boson Hamiltonian approach, in which Hit is 
taken as a collection of harmonic oscillators with various frequencies and the interaction 
term Hj is linear both in the j and in the ^ coordinates. One important result from this 
approach is that a two-level system & will describe damped oscillations between the localized 
states \R) and \L). Depending on the frequency distribution of the environment, & may 
be localized at T = 0°K (the overdamped case, known as "subohmic"), it may present 
critical damping (the "ohmic case") or it may undergo underdamped coherent oscillations 
(the "superohmic case.") The last one of these three instances is the most interesting for 
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the present discussion, as it allows the observation of MQC before the complete relaxation 
of the system. The possibihty of experimental MQC in the superohmic case depends in the 
interplay between a decoherence time defined only by the bath parameters, and the period 
of the Leggett-Caldeira oscillation for system (3. 

VI. THE SCALES OF MQC 



b 


Eq 


El 


AE 
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( nm ) 


(X 10-26 J ) 


( X 10-26 J ) 


(X 10-28 J ) 


( ) 


100.00000 
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2.9 


136.07900 
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5.4160961 


1.7 


3.8 


158.74011 


5.4036276 


5.4113353 


0.77 


8.6 


185.17494 


5.4059909 


5.4089897 


0.30 


22.0 


216.01195 


5.4069931 


5.4079902 


0.10 


66.0 


251.98421 


5.4073539 


5.4076298 


2.7x10-2 
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TABLE L Period r increases exponentially as the barrier width is augmented, a = 1.0 /xm, 
k = 2x 10-20J. This table, as well as all fi gures, was generated with Matlab ® R2012a. 

Let us start by fixing the width of the lateral wells at: 

a — l/im, (86) 

a value typical of contemporary lithographic circuitry, and take m to be the rest mass of an 
electron: 

m = me = 9.1 X IQ-^^kg. (87) 

With this, B takes the value: 

S = 0.6 X 10-2^J = 0.36 //eV, (88) 

and Tb is fixed at: 

Tb ^ 1.1 mK. (89) 



24 



From equation (l25|l . that gives the fundamental frequency of the Caldeira-Leggett oscil 
lations, we get the corresponding period 



El — Eq 

A global lower bound for this period is found from expressions f lT^ and f[T5]) : 



(90) 



2'Kh Ama^ , . 

For values (!86|) and (IHTl) this gives 

r > 11ns. (92) 

From table I (obtained through computer assisted numerical analysis) we get that as we 
sweep the barrier width from 0.2 to 0.5 yum the period of the Leggett -Caldeira oscillations 
for our square double well increases from 1.0 to 240 /is. Based on general considerations it 
has been estimated^ that, for all practical purposes, MQC is lost if the period of the Legget- 
Caldiera oscillation is of the order r > 100/is. Thus, the last row of the table corresponds to 
a localized system. All the other tabulated values could in principle correspond to observable 
MQC. 



A. Some of the many things we have left out 

MQC experiments are carried out in superconducting quantum interference devices 
(SQUlDs) with low capacitance tunneling Josephson junctions,-i^ and the relevant coordi- 
nate [i. e. the analogous of coordinate x) is not of a geometric character (like a position) 
but is in most cases the phase difference between the states functions of the electrons in a 
Cooper pair (so that m is not really the mass of the electron.) Thus our toy model is in 
reality a simplification of a mechanical analogy used to discuss experimental MQC. 



VII. CONCLUSIONS 



Contemporary quantum mechanics, both experimental and theoretical, provides exam- 
ples of basic concepts and techniques such as: tunneling, stationary states, two-level systems, 
perturbation theory, the density matrix and the WKB approximation. Classroom presen- 
tations of current areas of research, such as MQC, help to improve the understanding of 
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quantum physics at university level, as they connect the simplified textbook models with 
the actual state of the field, and thus with the future professional activity of the student. 
Moreover, MQC illustrates in a beautiful way the interplay between theory and experiment, 
and between concepts and techniques arising in different areas of quantum physics. 

We believe to have achieved in the present paper a level of exposition that makes it both 
clear and interesting for senior university students and recent graduates. In order to do 
so, we had to glide over the more technical aspects of experimental MQC and the intricate 
relation between MQC and the epistemology and the philosophy of physics. We hope that 
the present paper will encourage the interested reader to delve further into this facets of 
contemporary research. 



VIII. APPENDIX 



Consider condition (|8]) for the ground level (n = 0), that is: 

E, coe = {k- Eo) tanh^ ^ VMfc-i^o) (93) 



We will now establish a lower bound for Eq starting from (l93l) . but we have to take some 
precautions in doing so because Eq depends implictly on b. In order to proceed, note that 

V6.(0.oo), VM|EM<v{M*EM (94) 

ft ft 

so that 



^,J2m{k-Eo) ^ ^ .J2m{k-B/A) 



V6 G (0, 00) , tanh^ ^- > tanh^ — (95) 

It/ ft 

The dependence of the rhs of inequality (!95|) is explicit, so that the usual procedures of 
calculus can be applied. In particular as we now from elementary theorems that the limit 

lim tanh^ ^y^2^(k^Bj^ ^ ^ 

fe-s>oo h 

holds true, we can affirm that: for given 6 > there exists a bo^S) such that any b > bo^S) 



tanh^ b^ ^— '-^ > 1 - ^ (97) 

h 2k ^ ' 

From fl93l) fl95l) and fl97|) we deduce that for any b above a certain value bQ{6), the ground 

energy of Ub satisfies: 



Eocot^a^^>ik-Eo){l~^) (98) 
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Turning our attention to the condition for Ei, i. e. 

cot^ = ik- E,) coth^ 5 (99) 

we now find an upper bound for Ei, by noting that, because of (194|) and the known properties 
of the hyperbohc functions, the inequahty 



coth^ ^- < coth^ ^— — (100) 

n h 



is verified for all strictly positive b. Furthermore, 

itocoth^6^MEEM = i (101) 

SO that for every 6 > there exist a bi{6) such that, if 6 > bi{6), then inequality 



^,^ v/2m(fc-i?/4) 5 



coth^ 6 V 1^ < 1 + -t: (102) 

is satisfied for all strictly positive b. And from fl99l) and fll02p we get that, for all b above a 
certain thershold value bi{6), the inequality 



E,coeay^<ik-E,){l + ^) (103) 

is satisfied. 

Taking both (l98l) and (11031) into consideration, we have that for every 5 > there exists 
a number b{6) = max{6o(5)&i(5)} such that for any b > b{6) the inequality 



E, cot' - Eo cot' < [E, - E,) ^ bi\ - ^^±ll\ (104) 

a n \ 2k / 



is satisfied. Now, it is not difficult to see that 



2 V2mE 



w{E) = E cot^ a^-^ — (105) 
is a monotonically increasing function of E in the range B/A < E < B, so that 



2 \i2imE[ 2 V2mEo 



0< El cot^ a ^ ^ - ^0 cot^ a ^ (106) 
ft ft 



and in the other hand, we deduce 

Eq + Ei\ ( Eq + El" 



,,„.,,,s{^-^)<{^-^y<s COT) 
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from and . From {MM, ffT06|) and ffTOTj) we get: 



2 ^j2mEi 2 „ V^mEo 



0< El cot^ ^ - Eo cot^ a " ^ ^ <S (108) 
n h 



El cot^ a^L?p^ > k-Ei (109) 

/ 1 



Finally, we notice that, as 



, V2mE 



v(E) = cot^a—- — (110) 

h 



is a monotonically increasing function of E in the range B/4 < E < B, then 

(Ei-Eo)cot2a^^^?|5<5 (111) 

Now we just need to find a lower bound on cot^ a^^^^^y^. This is obtained by turning back 
to condition (p3|) from which we get 

cot^a ^ ^ ° < '— 112 

« B 

Finally, from fillip and f lll2p we arrive at 

Let us stress that h and 5 are independent of h. In this manner, we have arrived at the 
following lemma: 

For each strictly positive real number 5 there exists a 

6'(5) = 6(5^^1^) (114) 

such that for any h > b'{6) the gap between the ground and first excited levels of of Ub is 
less than 6, that is, such that: 

Ei-Eo<S. (115) 
And this is what we set out to prove in this appendix. 
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